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We report ab-initio calculations of the phonon dispersion relations of the single-layer and bulk 
dichalcogenides M0S2 and WS2. We explore in detail the behavior of the Raman active modes, 
Aig and as a function of the number of layers. In agreement with recent Raman spectroscopy 
measurements [C. Lee et. al., ACS Nano 4, 2695 (2010)] we find that the Aig mode increases in 
frequency with increasing layer number while the mode decreases. We explain this decrease by an 
enhancement of the dielectric screening of the long-range Coulomb interaction between the effective 
charges with growing number of layers. This decrease in the long-range part over-compensates the 
increase of the short-range interaction due to the weak inter-layer interaction. 
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I. INTRODUCTION 

Recently, the mechanical exfoliation technique, ap- 
plied to layered materials, has lead to the fabrication of 
new bidimensional systems, i.e., atomically thin layers.^^ 
Graphene,^ a planar sheet of carbon atoms arranged in 
a hexagonal lattice, is the most famous bidimensional 
material and exhibits intriguing physical properties not 
found in its bulk counterpart graphite.'^ However, the 
absence of a bandgap makes its use in electronic devices 
(transistors) difficult. Several strategies were proposed 
to overcome this setback by opening a gap: quantum 
confinement in nanoribbons,^ deposition of a graphene 
monolayer on boron nitride,^ applying an electric field 
in bilayer graphene,^ etc. Nevertheless, the experimental 
realization of a bandgap larger than 400 meV remains a 
challenge,!^ besides of a deterioration of other graphene 
properties, in particular the high mobility. Therefore, the 
fabrication of atomically thin sheets of other materials, 
with a finite bandgap, appears as the natural strategy in 
the search of materials for a new generation of electronic 
devices. 

In recent experiments, molybdenum disulfide (M0S2), 
an indirect semiconductor of bandgap 1.29 eV in its bulk 
phase, has exhibited a direct bandgap of 1.75 eV in 
its single-layer phase, together with an enhancement of 
the luminescence quantum yield in comparison with the 
M0S2 bulk.^ Moreover, Radisavljevic et al^ have 
demonstrated suitable properties of a single-layer M0S2- 
based transistor, like a room-temperature electron mo- 
bility close to that of graphene nanoribbons and a high 
on/off ratio. Therefore, single-layer M0S2 has become an 
appealing material in the area of optoelectronic devices, 
being an alternative and/or complement to graphene. 
From other layered compounds such as WS2, M0S2, 
BN,. . . monolayers can be produced by (liquid) exfolia- 
tion as welLl^ Moreover, M0S2 and other layered mate- 
rials are also interesting due to change of certain proper- 
ties with respect to its bulk counterparts. Finally, their 
topology facilitates the chemical identification atom-by- 
atom. 1^ 



Recent Raman spectroscopy measurements of M0S2 
single and multi-layers have revealed unexpected trends 
of the vibrational properties when the number of layers 
changes. Lee et. ali^ reported a decrease in frequency 
of the optical E^g phonon mode with increasing num- 
ber of layers. This is consistent with the finding that in 
bulk M0S2 the infrared active Eiu mode (where neighbor- 
ing layers are vibrating in-phase) is slightly lower in fre- 
quency than the Raman active E^g mode (where neigh- 
boring layers are vibrating with a phase shift of ir)^ 
But both observations contradict the naive expectation 
that the week inter-layer forces should increase the ef- 
fective restoring forces acting on atoms. One would thus 
rather expect a slight increase of the frequency of the Ra- 
man active mode with respect to the IR active mod^^^ 
and, accordingly, an increase of the frequency of the bulk 
Raman active mode with respect to the corresponding 
single-layer mode. As a plausible explanation of this 
anomalous trend the long range Coulomb interaction was 
mentioned.'^ The purpose of our article is the clarifi- 
cation of this issue by a detailed ab-initio study of the 
inter-atomic force constants, separating the short-range 
and the long range Coulomb parts. We show in the fol- 
lowing, that the anomalous trend in the £^2^ mode fre- 
quency is caused by the dielectric screening of the long 
range Coulomb forces in bulk M0S2. At the same time, 
we present a full ab-initio study of the phonon dispersion 
relations of the phonon dispersion relations of single-layer 
and of bulk M0S2 and the intimately related material 
WS2. (Tungsten is in the same column of the periodic 
system as Molybdenum). Apart from a study of the vi- 
brational stability of M0S2 nanoribbons,!^ a fully com- 
prehensive ab initio study of the vibrational properties 
of these materials is still absent in the literature. 

In Section II we present the methods for the calcula- 
tion of the phonon dispersions and the analysis of the 
force constants. In section III we discuss the phonon dis- 
persion relations of M0S2 and WS2 single layers and bulk 
and compare with experimental data. In section IV, we 
present the calculated results for the evolution of the Ra- 
man active phonon modes as a function of the number of 
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FIG. 1. (a) M0S2 bulk and single- layer. The interlayer dis- 
tance is denoted by d. (b) and (c) Slide and top view of the 
M0S2 bulk unit cell (analogously for WS2) unit cell. The 
primitive vectors are a = (a, 0,0), b = (a/2, v^a/2, 0), and 
c = (0,0,c). The weak layer interaction is indicated by a 
spring. 



layers and give an explanation in terms of the short-range 
and long-range contributions to the force constants. 



II. CALCULATION METHODS 

M0S2 and WS2 belong to the dichalcogenide family of 
materials, built up of weakly (van der Waals) bonded S- 
Mo-S single-layers as shown in Fig. [l] Each one of these 
single-layers consists of two hexagonal planes of S atoms 
and an intercalated hexagonal plane of Mo atoms bound 
with the sulfur atoms in a trigonal prismatic arrange- 
ment. The symmetry space group of bulk M0S2 is P3ml 
(point group Dq^)- The space group of the single layer 
is P6m2 (point group Dsh)- Consequently, systems with 
even number of layers belong to the space group P3ml 
(with inversion symmetry), and systems with odd num- 
ber of layers to the P6m2 space group (without inversion 
symmetry). 

The phonon calculations begins with the determina- 
tion of the equilibrium geometry (i.e. the relative atomic 
positions in the unit cell which yield zero forces and 
the lattice constants which leads to a zero-stress ten- 
sor). The calculations have been done with density 
functional theory (DFT) as implemented in the open- 
source code ABINIT,!^ within the local density approx- 
imation (LDA).^^We use Hartwigsen-Goedecker-Hutter 
pseudopotentials^^ (including the semi-core electrons as 
valence electrons in the case of Mo and W) and a plane- 
wave cutoff at 60 Ha. The first Brillouin zone is sampled 
with al2xl2x4 Monkhorst-Pack grid for bulk and 
12 X 12 X 1 for single and few- layers systems. 

The optimized lattice parameters are shown in ta- 
ble [l[ The experimental lattice parameters of M0S2 are 
a = 3.15 and c = 12.3 A.^ In the case of WS2 they are 



a = 3.153 and c = 12.323 A.'^^' Qur LDA calculations 
underestimate them by 0.7 % and 2.1 %, respectively. 
The slight underestimations of the in-plane lattice con- 
stant is a common feature of the LDA which tends to 
overestimate the strength of covalent bonds. For the 
single-layer, we checked the influence of the exchange- 
correlation potential on the geometry and the phonon 
dispersion by performing calculations within the Gener- 
alized Gradient Approximation (GGA).^For the single- 
layer of both M0S2 and WS2, we obtain a lattice constant 
a = 3.18 A, 1.76 % larger than the LDA value and 0.96 % 
larger than the experimental (bulk) value. Correspond- 
ingly, the phonon frequencies are reduced by an almost 
constant factor of 1.04 % throughout the whole phonon 
dispersion. 

The correct description of the c parameter is less ev- 
ident because the LDA (and other semi-local function- 
als) completely neglect the van der Waals component 
of the inter-layer-interaction. At the same time, how- 
ever, the LDA strongly overestimates the (weak) cova- 
lent part of the inter-layer bonding. Thus, the LDA 
has quite successfully reproduced the geometry and also 
given reasonable results for layer phonon modes of dif- 
ferent layered materials such as graphite^^ and hexago- 
nal boron nitride (hBN),^^ as well as of single-layers of 
hB]>P^ and graphene^^ on a Ni(lll) surface. We thus 
expect (and the obtained value for c supports this ex- 
pectation) that the LDA works reasonably well for the 
inter-layer phonons of the M0S2. We note that a physi- 
cally correct description of the equilibrium geometry, and 
the potential energy surface around, would require the 
proper treatment of van der Waals contribution, e.g. on 
the level of the random-phase approximation as it has 
been done for bulk hBN^^ and for graphite. Alterna- 
tively, non-local functionals that are optimized for the 
description of the van der Waals interact ioiP^l could be 
used. Both approaches are, however, out of the scope of 
the present work. 

For the calculations of single-layer and few-layers sys- 
tems, we have used a periodic supercell, leaving enough 
distance between adjacent sheets. For instance, we use 
c = 13.25 Angstroms in the case of a single-layer. The re- 
maining interlayer interaction has negligible effects on the 
phonon frequencies. All the results show a slight reduc- 
tion of the in-plane lattice constant together with a slight 
stretching of the vertical distances, with the total effect 
of a smaller Mo-S bond length for decreasing number of 
layers, being 2.382 A for single-layer and 2.384 A for bulk 
M0S2. 

Once the equilibrium geometry has been obtained, the 
phonon frequencies u can be calculated. These phonon 
frequencies are the solution of the secular equation 
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where q is the phonon wave- vector, and Mj and Mj are 
the atomic masses of atoms / and J. The dynamical 
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Lattice constants (A) 

M0S2 1-layer 2-layer Bulk 

a 3.125 3.126 3.127 

c 12.066 

WS2 l-layer 2-layer Bulk 

a 3.125 3.125 3.126 

_c 12.145 

TABLE L Equilibrium lattice parameters of M0S2 and WS2 
obtained in this work. 



matrix is then defined as 
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dur{q)du^j{q) 



(2) 



where uj{q) denotes the displacement of atom / in di- 
rection a. The second derivative of the energy in Eq. |2] 
corresponds to the change of the force acting on atom / 
in direction a with respect to a displacement of atom J 
in direction 
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The Fourier transform of the gr-dependent matrix 
leads to the real space atomic force constant matrix 
Cia,j(3{Rij)^ where Rjj is the vector that joins atoms 
/ and J. Thus, Cia,j^ < (> 0) means a binding 
(anti-binding) force in direction a acting on atom / when 
atom J is displaced in direction p. It is worth to mention 
that the diagonal term in the atom index, Cia,i(3^ cor- 
responds, according to Newton's third law, to the total 
force exerted in the a-direction on the atom /, when the 
displacements of the atoms J in the /3-direction is are set 
to oneP^ 
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This term is always positive (unless the crystal is 
unstable) and in the following we refer to it as self- 
interaction. Eq. Q demonstrates the contribution of 
many atoms to the self-interaction. One can distinguish 
two contributions, the short range part (which is mainly 
due to covalent bonding to the close neighbors) and the 
Ewald or long range part^^ (due to the Coulomb forces 
between the effective charges). This distinction will be 
helpful to interpret the evolution of the self interaction 
for varying layer thickness and to understand the unex- 
pected trends of the phonon frequencies (section [lV|) . 

For the calculation of the dynamical matrix we have 
used density functional perturbation theory (DFPT)^^ 
where atomic displacements are taken as a perturbation 



potential and the resulting changes in electron density 
and energy are calculated self-consistently through a sys- 
tem of Kohn-Sham like equations. Within this approach 
the phonon frequency can be obtained for arbitrary gr, 
with calculation only in a single unit-cell. 

Since M0S2 and WS2 are slightly polar materials, cer- 
tain IR active phonon modes at F give rise to a macro- 
scopic electric field. This electric fields affects the longi- 
tudinal optical (LO) phonons in the limit g ^ 0, break- 
ing the degeneracy of the LO mode with the transversal 
optical (TO) mode.lSSlThus, in bulk M0S2 and W0S2, the 
nonanalytic part of the dynamical matrix (which contains 
the effective charges and the dielectric tensor) must be 
calculated in order to obtain the correct frequencies at 
the Brillouin zone-center.^^ The LO-TO splitting for the 
Eiu mode has the value of 2.8 cm~^. In the case of a 
single or few-layers system, this effect is even smaller. 



III. PHONON DISPERSIONS 
A. M0S2 

We start our analysis of the vibrational properties with 
the description of the general features of the phonon dis- 
persions of bulk and single-layer M0S2, shown in Fig- 
ure |2] We have also depicted the experimental data ob- 
tained with neutron inelastic scattering spectroscopics. 
The overall agreement between theory and experiment is 
good, even for the inter-layer modes. This confirms our 
expectation that the LDA describes reasonably well the 
inter-layer interaction (even though not describing the 
proper physics of the inter-layer forces). 

The bulk phonon dispersion has three acoustic modes. 
Those that vibrate in-plane (longitudinal acoustic, LA, 
and transverse acoustic, TA) have a linear dispersion and 
higher energy than the out-of-plane acoustic (ZA) mode. 
The latter displays a g^-dependence analogously to that 
of the ZA mode in graphene (which is a consequence of 
the point-group symmetr}B21) ^ xhe low frequency optical 
modes are found at 35.2 and 57.7 cm~^ and correspond 
to rigid-layer shear/ vertical motion, respectively (in anal- 
ogy with the low frequency optical modes in graphite^Sl) . 
When the wave number q increases, the acoustic and low 
frequency optical branches almost match. It is worth to 
mention the absence of degeneracies at the high symme- 
try points M and K and the two crossings of the LA and 
TA branches just before and after the M point. 

The high frequency optical modes are separated from 
the low frequency modes by a gap of 49 cm~^. We have 
drawn in Fig. [3] the atomic displacements of the Raman 
active modes (£^2c/ ^1^) infrared active mode 

Eiu. The Raman active modes are also indicated in the 
phonon dispersion of Fig. |2j The in-plane modes 
and Eiu are slightly split in energy (by 3 cm~^). This is 
known as Davydov splitting and, for M0S2, the experi- 
mental value is 1 cm~^.^^ However, the finding that the 
Eiu frequency is larger than that of the E^g mode con- 
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FIG. 2. Phonon dispersion curves and density of states of 
single-layer and bulk M0S2. Points are experimental data ex- 
tracted from Ref. 18 Low panel, inset of the phonon branches 
in the region of the and Aig modes. 



tradicts a priori what one would expect from the weak 
inter-layer interaction: As one can see in Fig. [3j for the 
mode £^2^, the sulfur atoms of different layers move in 
opposite direction and thus the additional "spring" be- 
tween sulfur atoms of neighboring layers should increase 
the frequency of the mode £^2^ with respect of that of 
Eiu mode where sulfur atoms of neighboring sheets are 
moving in phase and thus the additional "spring" has 
no effect. The semi-empirical model of Ref. 16 takes 
this consideration into account, and obtained indeed that 
uj^i > ujeiu while experiment^^^ demonstrate the op- 
posite behavior. Our ab-initio calculations match the ex- 
perimental results which indicates that other causes be- 
yond the weak inter-layer interaction are present in the 
system. We will analyze this feature in the next Section 
with the aid of the atomic force constants. 

We now turn to analyze the single-layer phonon dis- 
persion, shown in Fig. |2j The symmetry is reduced from 
to Dsh and there is no longer a center of inversion as 
in the bulk. Therefore, the phonon mode labels at F must 
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FIG. 3. Phonon modes in-plane £2^, Eiu, and the out-of- 
plane phonon mode Aig, for the bulk Mo5'2 (analogously for 
WS2). 
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TABLE 11. Relevant phonon symmetry representations of 
single- layer (point group Dsn) and bulk (point group Dqh) 
M0S2 (inspired in Table II of Ref. [T5|) . Direction out-of-plane 
(in-plane) is parallel (perpendicular) to the c vector of the 
unit cell, respectively. Phonon frequencies are the calculated 
values of this work. 



be changed accordingly. The number of phonon branches 
is reduced to nine. Table HIl shows the most relevant 
M0S2 single-layer and bulk modes at F, together with 
their character, displacement direction, involved atoms, 
and frequency. 

Overall, the single-layer and bulk phonon dispersions 
have a remarkable resemblance. In the bulk, all single- 
layer modes are split into two branches but since the 
inter-layer interaction is weak, the splitting is very low 
(similar to the situation in graphite and graphene.'^S The 
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only notable exception from this is the splitting of the 
acoustic modes around P. In the single-layer, the result- 
ing low frequency optical modes are absent. 

In the single layer, the high frequency T modes 
and Eiu collapse into the mode E' . (From Fig. [s] it 
is evident that with increasing inter-layer distance, the 
modes £^2^ and Eiu acquire the same frequency.) Inter- 
estingly, as measured in Ref. 14 and indicated in Table [TH 
the bulk E2g mode is lower in frequency than the single- 
layer E' mode. This contradicts the expectation that 
the additional inter-layer interaction should increase the 
frequency but is in line with the anomalous sign of the 
Davydov splitting between the bulk £^2^ and Eiu modes. 
The origin of this will be discussed in section |IV[ The 
out-of-plane mode Aig follows the expected trend that 
the inter-layer interaction increases the frequency with 
respect to the single-layer mode Ai. 

The densities of states (DOS) of single-layer and bulk 
are represented in the right panels of Fig. [3] The differ- 
ences between single-layer and bulk DOS are minimal, 
except a little shoulder around 60 cm~^ in the bulk DOS 
due to the low frequency optical modes. In both cases 
the highest peaks are located close to the Raman active 
modes E2g and Aig. 



B. WS2 

Figure [4] shows the phonon dispersions of single-layer 
and bulk WS2, together with the density of states (DOS). 
The general features are identical to those of the disper- 
sions of M0S2 (Fig. [2|. The differences between single- 
layer and bulk dispersions are similarly weak as in the 
case of M0S2. Thus, also the bulk DOS resembles very 
much that of the single-layer (except for the small shoul- 
der ~ 50 cm~^ due to the inter-layer optical modes). 

For a better comparison of M0S2 and WS2 single-layer 
phonon frequencies, we have depicted them together in 
Fig. [5J In general, the WS2 phonon bands are shifted 
down to lower frequencies with respect to the M0S2 fre- 
quencies. The cause of this trend is the larger mass of 
the tungsten atoms, and therefore their lower vibration 
frequency (see Eq. [T]) . The only notable exceptions from 
this general downshift are the branches associated to the 
mode E''^ around 300 cm~^ and to the Ai mode around 
410 cm~^. In these modes, only the sulfur atoms are 
vibrating (see Table |ll| and thus their frequency is not 
affected by the mass of the metal atom (W or Mo), just 
by the strength of the covalent bond. 

The larger mass of W leads to a larger frequency gap 
between low and high frequency modes (110 cm~^) since 
the highest acoustic branch is pushed down. Further- 
more, the difference between the modes Ai and E' ^ is 
now of 60 cm~^, three times larger than in the case of 
M0S2. 

The density of states of the WS2 single-layer is also 
appreciably different from that of M0S2. While at low 
frequencies the DOS has two well differenced peaks, as 
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FIG. 4. Phonon dispersion curves and density of states of 
1-layer and bulk WS2. 




FIG. 5. Phonon dispersion curves for single- layer of M0S2 
(solid lines) and WS2 (dashed lines). The density of states is 
depicted in the right panel. 



in Fig. [2j for higher energies one peak stands out from 
the others, at frequency of ~ 350 cm~^, and associated 
mainly to the P-point mode E' . 



IV. EVOLUTION OF Aig AND E^g PHONON 
MODES WITH THE NUMBER OF LAYERS 



The understanding of the frequency trends of the Aig 
and E2g modes with varying layer thickness requires a 
more refined analysis. With the aim of explaining the 
Raman scattering experiments of Ref. T!4^ we have calcu- 
lated the phonon frequencies at the P point for single, 
double- and triple-layers and we discuss the evolution of 
the atomic force constants from single-layer to bulk in 
detail. 

Figure [6] shows the phonon frequency of Aig and E^g 
modes as a function of the number of layers. Since LDA 
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1 2 3 4 5 6 Bulk 
n-layers 



FIG. 6. Phonon frequencies of Aig and modes as a 
function of M0S2 layer thickness, as obtained in this work 
(squares) and experimental data from Ref . ^2 (circles) . We 
plot the frequency differences with respect to the correspond- 
ing bulk modes. The insets represent the phonon modes at F 
point. 



tends to overestimate the phonon frequencies, it is rea- 
sonable to represent the difference between the n-layers 
frequency and the bulk frequency instead of comparing 
absolute theoretical and experimental values. In such a 
representation we observe that our calculation properly 
reproduce the up-shift of the Aig mode and the down- 
shift of the mode with increasing layer number. The 
quantitative differences between theory and experiment 
are mainly due to the limited precision of the description 
of the inter-layer interaction by the LDA. 

After we have shown that DFPT-LDA reproduces the 
experimental trends we would like to give an explana- 
tion of this behavior, taking advantage of the detailed 
knowledge of the atomic force constants, available in ab 
initio methods. We interpret the changes in phonon fre- 



quency through an analysis of the real-space force con- 
stants (Eq. |3|, in particular of the self-interaction term 
(Eq.[4|. Going from the single-layer to a multi-layer, this 
term changes in two ways: i) the short range term in- 
creases due to the (weak, but non-negligible) interaction 
with atoms from neighboring layers; (ii) the long-range 
Coulomb interaction changes, because the sum extends 
over effective charges from all layers and the effective 
screening of the Coulomb potential increases. 

We start our analysis with the out-of-plane mode Aig. 
From Fig. [3j it is intuitively clear that the interaction 
between sulfur atoms of neighboring layers can influence 
substantially the frequency. Going from the single-layer 
to the double-layer, one "adds" an additional "spring" 
between the atoms S and S' on neighboring layers which 
leads to an increase of the Aig mode frequency with in- 
creasing number of layers. As only sulfur atoms are in- 
volved in this mode, we just need to examine the sulfur 
self-interaction term {Csz,Sz) and the atomic force con- 
stant {Csz,S'z) between nearest neighbors, that belong 
to adjacent layers (atoms joined by the spring in Fig. [3|. 
We have summarized in Table [TTTl the results. Note that 
the small variation of the term Csz,Sz from single-layer to 
bulk prevents it from being the main cause of a frequency 
increasing of almost 5 cm~^. Thus, the term Csz,S'z has 
a negative value which implies a binding force ( "spring" ) 
between atoms S and S' which leads to an increase of 
frequency. Force constants related to farer neighbors are 
negligible in comparison with Csz^S'z- This demonstrates 
that the weak interlayer interaction is the main cause of 
the frequency increasing with the number of layers. 

One might expect that the same argument holds for 
the E2g mode: the additional "spring" between sulfur 
atoms from neighboring layers should increase the fre- 
quency with increasing number of layers as well. How- 
ever, theoretical and experimental results show the op- 
posite behavior. The reason can be found in the self- 
interaction terms CmoxMox 

for Mo and S in Table [Sj 
while the difference between single-layer and bulk is neg- 
ligible for the sulfur atoms, one observes a considerable 
decrease for the molybdenum atoms. Interestingly, the 
short range contribution to the self-interaction increases 
(as one would expect from the inter-layer interaction). 
However, the long-range Coulomb part^^ decreases con- 
siderably. In the appendix, we show that this decrease 
is related to the strong increase of the dielectric tensor 
(both in-plane and out-of-plane) when going from the 
single-layer (represented in our calculations by a periodic 
stacking of single-layers with large vacuum between) to 
the bulk (see Table [TTT| ) . We note that one might asso- 
ciate the change in frequency to the differences in lattice 
constant and interatomic distances in bulk and single- 
layer, respectively. However, the small differences shown 
in Section II are not enough to account for the observed 
magnitude of the E^g frequency change. 

Therefore, the decrease of the E^g phonon frequency 
is associated to a stronger dielectric screening of the long 
range Coulomb interaction in few-layer and bulk M0S2. 
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The effect is particularly pronounced for the molybde- 
num atoms, as Table |III| shows. Our analysis also ex- 
plains why previous empirical models have not been able 
to explain the experimental observations, due to the dif- 
ficulty to include this subtle change in the parameters. 



V. CONCLUSIONS 

In conclusion, we have studied the phonon dispersions 
of M0S2 and WS2 single layers and bulk using density 
functional perturbation theory in the local density ap- 
proximation. We obtain good agreement with neutron 
diffraction data as well as Raman and infrared absorption 
spectroscopies. We have explored how the Raman active 
modes Aig and E^g change their frequencies when the 
number of layers varies, and confirm the recently reported 
down-shift of the mode with increasing number of 
layers. This unexpected behavior can be explained by an 
increase of the dielectric screening which reduces the long 
range Coulomb interaction between the effective charges 
and thus reduces the overall restoring force on the atoms. 
We expect that this explanation also holds for other polar 
layered materials where an anomalous Davydov splitting 
has been observed (such as GaSe^^ and GaS^^) or where 
the Raman peak shifts down with increasing number of 
layers such as it has been recently observed for hexagonal 

BN.Ea 



and the dielectric tensor. The long range contribution, 
j^, can be written in terms of the dielectric tensor 
e, its inverse, €~^, the Born effective charges Zj and 
the interatomic distance d = djj^ as defined in Ref. |33l 

^Ir _ r^* r^* f {^~^)a'f3' ^ ^a'^P' \ 

where = X]/3(^~"'^)a/3^/3 conjugate of the vec- 

tor d, and the norm of the latter in this metrics is 
D = VA • d. This expression simplifies enormously by 
assuming diagonal the dielectric tensor, and Zj = 
Zj^^,Saa'' After some algebra one obtains: 

(Al) 

We can examine with an example how the long-range 
atomic force constants of single-layer and bulk M0S2 are 
related with the dielectric tensors. Thus, we can evaluate 
the term for neighbor Mo atoms that belong to 

the same layer, with interatomic distance d = (d, 0,0), 
and assuming the same distance for single-layer and bulk. 
Thus, we obtain a simplified expression of Eq. |A1[ 
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^Mo,x,Mo,x 



(A2) 



The Born effective charges Z]^^ are almost equal 
for both systems, and using the dielectric tensors given 
in Table IIIIl we obtain: 



Appendix A 

The examination of the long-range atomic force con- 
stant formula can help to establish quantitatively its 
change with the variation of the dielectric tensor in the 
single-layer and bulk M0S2. However, the anisotropy 
of the crystalline structure of M0S2 impedes a direct 
relation between the long range atomic force constants 



^Mo,x,Mo,x{^^) 
^Mo,x,Mo,xi^'^^^) 



^xx,bulk^zz,bulk 
^xx,ll^zz,ll 



3.09. (A3) 



From the ab-initio calculation of the atomic con- 
stants we obtain C]^o,x,Mo,xi^^)/^Mo,x,Mo,xi^^lf^) = 
3.19, which is in agreement with the value obtained in 
Eq. A3 This demonstrates that the difference in the 
long-range part of the force constants for single-layer and 
bulk originates from the different dielectric screening. 
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